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We demonstrate the design of an integrated conditional quantum teleporta- 
tion circuit and a readout circuit using a two-dimensional photonic crystal 
single chip. Fabrication and testing of the proposed quantum circuit can be 
accomplished with current or near future semiconductor process technology 
and experimental techniques. The readout part of our device, which has 
potential for independent use as an atomic interferometer, can also be used on 
its own or integrated with other compatible optical circuits to achieve atomic 
state detection. Further improvement of the device in terms of compactness 
and robustness could be achieved by integrating it with sources and detectors 
in the optical regime. 
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1. Introduction 

It has been shown theoretically by Bennett 1 , et al. that quantum teleportation of an un- 
known state can be achieved using Bell states. Since then several theoretical and experimen- 
tal schemes have been proposed, which include implementation in optical, NMR, and cavity 
QED systems 2-5 . Inspired by the recent paper of S.-B. Zheng 2 , in this letter we propose a 
photonic crystal (PC) based cavity QED implementation of the scheme described therein, 
as well as the integration of the device with an atomic readout circuit to detect the tele- 
ported state. Since our proposed device functions as a very simple conditional teleportation 
chip, where the success of the protocol is conditioned on a specific measurement result con- 
trary to unconditional one, it could be fabricated and tested with current or near future 
semiconductor processing technology and experimental techniques, respectively. 

An illustration of the photonic crystal chip that we propose is shown in Fig. 1, where a 
single mode cavity is integrated with two parallel waveguides. The actual size of the chip, 
which is designed in this letter, is larger than that of the illustration. The actual height (i.e., 
vertical) is 17 ^ a and the width (i.e., horizontal) is 79a, where a is the lattice constant. x\, X2, 
£3, £4 and y are 9a, 43a, 18a, 9a, and 4\/3a, respectively. Because of its large size, we design 
the chip in a 2D lattice to reduce the high computational cost of 3D crystal analysis. This 
is, however, a good approximation due to the ability of 2D photonic crystal defect modes to 
emulate those of 3D photonic crystals. This emulation issue is considered in detail in Ref. 
[6] and its advantages for photonic crystal based quantum logic gate design are discussed in 
Ref. [7], where we have used exactly the same lattice parameters. Thus, our 2D PC structure 
is chosen to emulate its 3D counterpart and is based on a triangular lattice of silicon rods 
with dielectric constant of 12 and lattice constant of 2.202mm. The radii of the rods and the 
defects are 0.175a and 0.071a, respectively. The trajectories of the atoms through the PC 
are depicted by the red and blue dashed lines in Fig. 1. 

Incorporated with detectors, the cavity part of the chip operates as a teleportation circuit 
and the parallel waveguides act as a readout circuit to detect the teleported atomic state 
(see Fig. 1). Below we use the terms, teleportation circuit (TC) and readout circuit (RC), 
respectively to refer those two integrated photonic crystal circuits. We choose the width of 
the former to be 43a and the latter to be 36a in order to avoid unwanted coupling between 
the two circuits. In the first part of this letter we describe the teleportation mechanism and 
in the second part we explain how to detect the atomic state 8-11 within our designed PC 
quantum circuit (PCQC). 

2. Teleportation Circuit 

The TC involves the interaction of two Rb atoms with a resonant high-Q (assumed to be 
larger than 10 8 to prevent cavity decoherence) single mode cavity and the detection of one 



2 



of the atoms using a detector. The total Hamiltonian (Jaynes-Cummings Model) for the TC 
under the dipole and rotating wave approximations is given by 7 : 

Hatom-cavity = ^ cr J z + Tiua ] a + % ^ |#(rj) I {a ] oL + acr+) (1) 

3 j 

where the summation is carried over two atoms, A and B, u> is the resonant frequency, o z 
is the ^-component of the Pauli spin operator and a± are the atomic raising and lowering 
operators, a and cv are the annihilation and creation operators, respectively. g(vj) is the 
atom-cavity coupling parameter at the position of atom j and can be written as 

g(r j )=g E(r J )/\E(r rn )\. (2) 

go denotes the vacuum Rabi frequency and r m is the position in the dielectric where the 
electric field energy density is maximum. Note that we implicitly assume in Eq. (2) that the 
atomic dipole moments, fi\ , of the atoms are aligned with the photon polarization. 

In the first stage of the conditional teleportation protocol, atom B (whose trajectory is 
depicted by the red dashed line in Fig. 1), initially in the excited state |1)b, traverses the 
cavity, which is initially in the vacuum state |0), with velocity Vb- Since at this stage atom 
A, whose state is to be teleported, is not in play, index j only takes value B in Eq. (1). The 
composite state for atom B and the cavity at time t then can be written as 

|tf (t)> = cosG B (t)|l) B |0) - ?sinG B (t)|0) B |l), (3) 

where 

Gj(t)= f\g(T S )\dT. (4) 

J to 

The cavity we have designed in this letter supports a transverse magnetic (TM) polarized 
monopole mode as illustrated in Fig. 1. Based on the block-iterative plane wave expansion 
method 12 , with 32 grid points per a, we find the normalized frequency of the cavity mode to 
be 0.3733c/a. By setting a = 2.202mm, we tune the resonant wavelength to 5.9mm. At this 
wavelength, the atomic dipole moment for the Rb atom is taken to be 2 x 10 _26 Cm. 7 

Using these parameters, we obtain Gs{t) — 9n/4 at time t\ = 51.6/is at 18a from the left 
edge of the TC, if Vb is set to 767.7m/s, and hence atom-cavity entanglement in Eq. (3). As 
atom B travels across the cavity, the time-dependent atom-cavity coupling parameter and 
the evolution of the probability amplitudes of states |l)s|0) (blue) and — i|0) B |l) (pink) in 
Eq. (3) are displayed in Fig. 2. In our calculations we neglect the insignificant effect of the 
mode tail on the probability amplitudes of the entangled state — less than 1% beyond 18a. 
Another even less significant source of error in our design is the coupling of the waveguides 
in the RC to the cavity 13,14 , which is simply circumvented by keeping X2 in Fig. 1 sufficiently 
large (i.e., 43a). 
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Having entangled atom B and the cavity, the second stage of the conditional teleportation 
protocol is to inject atom A (whose trajectory is illustrated by the blue dashed line in Fig. 
1), which is initially (i.e., at time t\) in an arbitrary (or an unknown) quantum state, 

\<p) A = c \Q) A + c 1 \l) A . (5) 
Thus at time t\, the state of the whole system is written as 

l¥>(*i)> = 4| (co|0)a + Cl |l) A )(|l) B |0) - t|0> fl |l». (6) 

Once atom A arrives at the detector, right before the detection, say at time t\ + 1 2 , the 
state of the whole has evolved into 

Mh + h)) = -L{ Co |0) A |l) B |0) - ic Q \0) B [cosG A (t 2 )\0) A \l) - ismG A (t 2 )\l) A \0)} 

+c 1 |l> fl [cosG A (t 2 )|l) A |0> - ismG A (t 2 )\0) A \l)} 

-i Cl \0) B [cosV2G A (t 2 )\l) A \l) - ismV2G A (t 2 )\0) A \2)}}, (7) 

which would reduce to Eq. (6) of Ref. (2), if |g(r,,)| were spatially uniform in Eq. (1). 

In the derivation of Eq. (7), the contribution of |g(rs)| in Eq. (1) is neglected, since atom 
B is sufficiently far from the cavity when atom A is injected. At times close to t%, although 
|g(re)| « |(yf(r A )|, neither atom alone nor combined can start the Rabi oscillation. 7 On the 
other hand, as atom A proceeds its coupling to the cavity mode increases, while the coupling 
of atom B decreases. This means that the effect of |g(rs)| becomes even less significant for 
times greater than t\. 

The third stage of the conditional teleportation protocol is the measurement of the state 
of atom A. If the result is |1),4; the combined state of atom B and the cavity becomes 

*(ti+t 2 ) 

The denominator in Eq. (8) is simply the norm of the unnormalized state in the numerator. 
We obtain G A (t 2 ) = 7tt/4 in Eq. (8) if we set V A = 987m/s, and thus 

\4>)b = c \0) B + Ci|l) B . (9) 

That is, the initial arbitrary (or unknown) state of atom A has been teleported to atom B, 
conditional on detecting atom A in state 1 1)^, which occurs with probability 1/4. The time- 
dependent coupling parameter for atom A and the evolution of the probability amplitudes 
of states — colO^lO) (shown in pink), ci|l)#|0) (shown in blue), and — «ci|0)b|1) (shown in 
yellow) for the unnormalized state (i.e., the numerator) in Eq. (8) are shown in Fig. 3. Note 
that as time elapses Eq. (8) transforms into Eq. (9). Note also that atom B arrives at 31a 
at time t\ + t 2 (i.e., 88.9/is). After an additional distance of 12a, the fourth stage of our 
conditional teleportation and detection protocol begins as atom B enters the RC at time 
^1+^2 + £3 (i-e., 123. 3«s). Next we study how to detect its state. 



-cosinG A (i2)|0>B|0> + c 1 cosG yi (t 2 )|l) B |0) - 
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3. Readout Circuit 

We can write the atom-maser Hamiltonian, under the dipole and rotating wave approxima- 
tions, as 15 

= ^ + + V + ) (10) 

where 

= 2/z 10 |E(r)| 
ft 

is the Rabi frequency and u m is the frequency of the maser, whose polarization is also 
matched with the atomic dipole moment, /iio, of the atom. In the following we will first 
assume a constant fi(r) = Qq to simplify the explanation of the detection mechanism, and 
then we will modify the analysis to apply to our proposed PCQC, where Q(r) is not constant. 

Assume that atom B in state \4>)b [see Eq. (9)] enters a uniform field region at time t 
and interacts during a time interval t±. Then its probability amplitudes at time t + £4 can be 
computed to be 16 

co(t + U) = {[i cos 6 sin(^) + cos(^)]c + i sin 6 sm^K^y^ 2 , (12) 

Zi Z Zi 

Cl (t + t A ) = {1 sin 6 sin(^)e-^*c + H cos 6 sin(^) + C os(^)] Cl } e -^ t4/2 . (13) 
where cos0, sin0, and A are defined as follows. 



A=y/{u-«) m )2 + nl (14) 
(uj- UJ m ) 

cosB = ; , 15 

A 

sin6 = -^. (16) 

After this interaction, atom B interacts with a second uniform field region for the same 
time interval £4 starting at time t + t^. The probability amplitudes, co(t + 2t4) and ci(t + 2£i), 
in this case can be also found from equations (12) and (13) by replacing c , c\, and t with 
co(t + £4), ci(t + t^), and t + £4, respectively. 

Note that one can also treat these two uniform field regions as a single uniform field region 
with an interaction time of 2£ 4 . In this letter, however, we also intend to show that our RC 
has the potential to be an atomic interferometer, like a Ramsey interferometer. 16 In order 
to hint how this could work, we design two separate oscillating field regions based on two 
parallel waveguides as discussed below. Implementation of a fully working interferometer, 
however, requires further consideration. 

After atom B has interacted with the two uniform field regions, we can implicitly write 
the final state of atom B at time t + 2^4, before detection as 

100 + 2t 4 ))s = (cqCqo + ciCio)|0) B + (c c i + c lCll )\l) B (17) 
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where is the probability amplitude of the atom being in state \j) at time t + t 4 if it is 
initially prepared in state For example if the atom is initially in state |1) then it evolves 
into Cio|0) + Cn 1 1) at time t + t 4 . 

Having interacted with both waveguides, atom B is finally detected by an ionization 
detector, as shown in Fig. 1. The probability of finding it in the excited state is 

Pi = Icol 2 !^! 2 + | Cl | 2 l Cll | 2 + [^c^ + c«c«]2(c 01 ct 1 ) W + [4 r) c« + c i) C ^ ) ]2(c 01 ct 1 )«. (18) 

Superscripts (r) and (i) represent the real and imaginary part, respectively. 
We can simplify Eq. (18) by writing Eq. (9) in the form 

\<f>) B = cos(#/2)|0> B + sin(0/2)e^|l) B (19) 

where we have implicitly set Cq = cos(#/2), cjj = 0, c(' = sin(6'/2) cos<^, and = 
sin(6'/2) simp. 

Then, at four different A = uj m — uj values we evaluate |c i| 2 , |cn| 2 , 2{cqic\ 1 ) <j '\ and 
2(coiC^ 1 )^- ) directly and measure Pi. Thus having obtained four equations, we can find |co| 2 , 
|ci| 2 , [c^cf^ + Cq a? ], [co c i ~ c o c i ] an d hence \4>)b- Note that, although we have only 
two unknowns in Eq. (19), it is convenient to use four equations due to the form of Eq. (18). 
This demonstrates that one could detect the teleported state with our proposed device. In 
order to measure Pi, we need an atomic beam of identically prepared single atoms, since 
multiple observations are necessary to estimate the probability. Since the TC only has a 25 
percent conditional success rate, incoming atoms to the RC will not be identical. Thus if 
atom A of the relevant pair is not detected in the excited state, the corresponding atom B 
measurement is discarded from the calculation of Pi. 

Finally, we describe how to apply these ideas to our proposed PCQC. Remember that 
atom B in state \4>)b enters the RC at time t\ + 12 + £3. The steady-state mode profiles of 
the waveguides are calculated by two-dimensional-finite-difference-time-domain (2D-FDTD) 
method with discretizations of ^ and in the horizontal and vertical directions, re- 

spectively. We have designed the waveguides as coupled-cavity- waveguides 17 ' 18 to allow the 
atoms sufficiently large void regions through which to travel freely without a resonant-dipole- 
dipole-interaction or the Casimir- Polder effect (see Fig. 1). The steady-state cross section of 
the electric field magnitude, for A = 0, along the path of atom B is shown in Fig. 4. At the 
steady state the coupling between the two parallel waveguides 19-22 is not significant. 

In our calculations we assume that atom B interacts with two separate nonuniform waveg- 
uide modes of total width 36a, shown by the pink curve in Fig. 4. In other words we truncated 
their tails (shown by blue in Fig. 4.), which overlap with the cavity mode. The effect of this 
truncation in our design is observed to be negligibly small. We have observed that the cou- 
pling of waveguide modes on the cavity mode is less than the cavity coupling to waveguides, 
which does not effect the probability amplitudes more than 0.6%. 
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Since there is no feedback mechanism, as there is in the atom-cavity interaction, the atom- 
waveguide interaction can be described by a Markoff approximation. In other words, when 
an atom emits a photon into the waveguide, it leaves the interaction region immediately 
and cannot act back on the atom. Then the waveguides must be designed in such a way 
that total rate of spontaneous emission into the waveguide and the free-space lossy modes 
is sufficiently low as atom B traverses them. In our design considering estimated effective 
transverse cross section of corresponding three-dimensional waveguide mode the maximum 
probability of spontaneous emission for a given atom is estimated to be on the order of 
10~ 4 . 23 Thus, it is safe to assume that the Bloch vector of atom B can be manipulated 
without being correlated with the defect mode 24 of the waveguides. 

The total interaction time with the parallel waveguides is also 2t±, as in uniform case above, 
with t 4 (i.e., 51.6/is) for each. Note that, since the Rabi frequency, fi(r), is proportional to 
|E(r)|, the pink part of the curve in Fig. 4 also describes the position dependence of the Rabi 
frequency if multiplied by 2fi W /h. Thus the evolution of the teleported state, as it interacts 
with the parallel waveguides, can be determined numerically by exploiting equations (12) 
and (13) successively, t\ replaced with the maximum tolerated time-step, which is 44ns in 
our case (i.e. time for atom B to travel a/65). 

Thus, right before the detection, the final state is \§{t\ + t 2 + £3 + 2t 4 ))s, which can be 
written exactly in the form of Eq. (17). The probability of detecting atom B in the excited 
state is then given by Eq. (18), and following the same method as in the uniform case above, 
we can easily detect \<P)b, the state output by TC. 

As an example, assume that the final state of the TC [see Eq. (19)] is 



The corresponding calculated values of | C01 1 2 (green), |cn| 2 (red), 2(coic 41 ) < - r - 1 (light blue), 
2(coic 41 )**- ) (purple), and the value of P\ (dark blue), which will be measured, are shown 
in Fig. 5. Picking four different values of A, we can easily determine that 6 = ir/4 and 
<p = — 7r/6. 

4. Discussions 

Due to the electromagnetic scalability of photonic crystals one could also implement this 
device in more compact form in the optical regime. In this case, however, the electric field 
gradient and the momentum recoil of the atom would also have to be analyzed. In the 
microwave regime studied in this letter, the small momentum recoil due to a microwave 
photon does not cause external degrees of freedom (i.e., the position or momentum) of the 
atom to couple with its internal degrees of freedom (i.e., the ground or excited state). 

To implement our PCQC, one should also engineer the surface states, which reside at the 
air-photonic crystal interfaces. These surface modes may also become correlated with the 




(20) 
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atoms exiting the crystal. These modes can be handled, however, by appropriate termination 
of the crystal. 25-28 

5. Conclusions 

In conclusion, we have designed an integrated conditional teleportation circuit with readout 
in a single PC chip. Our proposed PCQC can be implemented by current or near future 
semiconductor processing technologies and experimental techniques. Even the detectors and 
sources illustrated in Fig. 1 could also be integrated — the feasibility of single atom detection 
on a PC chip has already been demonstrated 29 , as have PC based lasers. 30 ' 31 The RC part of 
our device can not only be used with the TC but also can be used independently, or integrated 
with other compatible optical circuits to read out atomic states. It also hints at the potential 
for implementation of atomic interferometers, such as the Ramsey interferometer, in photonic 
crystals. 

Achieving good mode-matching stability becomes quite cumbersome for many multiply 
nested interferometers in optical quantum computing. Photonic crystals, however, may be 
an especially promising paradigm as robust quantum circuit boards for delicate, next gener- 
ation, scalable optical quantum computing and networking technologies, as well as scalable 
quantum dynamic random access memory 32 , where qubits are refreshed using PC based 
cyclical quantum circuits. Integrated PC devices, such as planar waveguides and elementary 
interferometer modules, could replace the current bulky elements in linear optical quantum 
circuits. Furthermore, such devices could be made quite compact using PC integrated with 
single photon sources. 33 
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List of Figure Captions 

Fig. 1. Illustration of the proposed photonic crystal quantum circuit. This illustration is not 
to scale. 

Fig. 2. (a) The time-dependent atom-cavity coupling parameter and (b) the evolution of the 
probability amplitudes of states |l)s|0) (blue) and —^0)^11) (pink) in Eq. (3), as atom B 
travels across the cavity. 

Fig. 3. (a) Time-dependent coupling parameter for atom A and (b) the evolution of the 
probability amplitudes of states — co|0)b|0) (pink), ci|1)b|0) (blue), and — ici|0)s|l) (yellow) 
for the unnormalized state (i.e., numerator) in Eq. (8). 

Fig. 4. The steady-state cross section of the electric field magnitude, for A = 0, along the 
path of atom B (see text). 

Fig. 5. Calculated values of |c i| 2 (green), |cn| 2 (red), 2(c ic* 1 )W (light blue), 2(001^)^ 
(purple), and to be measured value of Pi (dark blue) (see text). 
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Fig. 1. Illustration of the proposed photonic crystal quantum circuit. This 
illustration is not to scale. guneyFl.eps. 
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Fig. 2. (a) The time-dependent atom-cavity coupling parameter and (b) the 
evolution of the probability amplitudes of states |l)s|0) (blue) and — i|0)s|l) 
(pink) in Eq. (3), as atom B travels across the cavity. guneyF2.eps. 
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Fig. 3. (a) Time- dependent coupling parameter for atom A and (b) the evolu- 
tion of the probability amplitudes of states — co|0)s|0) (pink), ci|1)b|0) (blue), 
and — zci|0)b|1) (yellow) for the unnormalized state (i.e., numerator) in Eq. 
(8). guneyF3.eps. 
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Electric Field Magnitude vs. Position 
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Fig. 4. The steady-state cross section of the electric field magnitude, for A = 0, 
along the path of atom B (see text). guneyF4.eps. 
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Fig. 5. Calculated values of | c i | 2 (green), |cn| 2 (red), 2(c ic* 1 p r > (light blue), 
^{coiCii)^ (purple), and to be measured value of Pi (dark blue) (see text). 
guneyF5.eps. 
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